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1. Introduction 

The spin-charge separation is a main subject in one-dimensional interacting electron 
systems. The conformal field theory has succeeded in the description of spin-charge 
separation in the low-energy physics of the Tomonaga-Luttinger liquid. Beyond the 
conformal field theory limit, exactly solvable models provide us chances to obtain 
analytical knowledge on thermodynamics and dynamics, and it is intriguing how the 
spin-charge separation appears in these properties. 

Among exactly solvable models, the supersymmetric t-J model with l/r^ 
interaction [J reveals the spin-charge separation in the simplest manner. The 
Hamiltonian of this model is given by 



i<j 



hY^s^, (1) 



where Ci^ = Cio-(l — ni^^„) with being the annihilation operator of an electron 
with spin a at site i, and ni = X^o- "^i-o- — X^o- '^L'^icr- The spin operator associated 
with site i is defined as Si = X)a /3 '^L('''i)ct/3Ci/3/2 where cr = {cTx, Uy, az) is the 
vector of Pauli matrices. The transfer energy and exchange one Jij are given by 
tij — Jij/'i, = tD^j^ where Dij — {N / n) sin {n {i — j) /N) with N being the number 
of lattice sites. Henceforth we take t as the unit of energy. We note that this t- 
J model reduces to the Haldane-Shastry spin chain model [3 El at half-filling. In 
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the supersymmetric i- J model with interaction, exact thermodynamics can be 
interpreted in terms of free spinons and holons At low temperature, the spin 

susceptibility is independent of the electron density n, and the charge susceptibility 
is independent of the magnetization m. These features are referred as to the strong 
spin-charge separation In addition, the fact that the magnetization m for a certain 
range of h is independent of n can be regarded as the strong spin-charge separation. 
Namely, m is determined by h as follows: 



1 - 1/1 -2/1/7^2^ for 0<h<hc 
a, for h> he 

where he = n{2 - n)7rV2 0. 

The strong spin-charge separation appears also in dynamics at zero temperature. 
The dynamical spin structure factor is given by 

S^\Q,iu) = ^ |(i.|5^|0)p5(a; - E, + Eo), (3) 

where Sq = Sf e^^^ / \/N . Here denotes a normalized eigenstate of the 
Hamiltonian with energy E^, (Eq being the ground state energy). In the absence 
of magnetic field (/i = 0), the dynamical spin structure factor was exactly obtained at 
n = 1 [7||S1|H]- It was numerically demonstrated that the weight of the dynamical spin 
structure factor in the t-J model does not depend on n in the region where only two 
spinons contribute |1(J| . This is an indication of the strong spin-charge separation in 
dynamics. Later, a mathematical poof was given to this statement, and the analytical 
expression of the dynamical spin structure factor for n < 1 was obtained in the full 
{Q,u}) space [TT] . 

A numerical study ^] also showed that the strong spin-charge separation for 
S^^{Q,uj) can be extended to the case of finite magnetic field {h 7^ 0). Namely, at 
fixed magnetization, S^^{Q, ui) away from half-filling is the same as that for half-filling 
(i.e., the Haldane-Shastry model), in the region where only spinons and antispinons 
contribute. For h 0, the full exact results on S^^{Q, cu) have not been obtained even 
in the Haldane-Shastry model. However, if the momentum is restricted to Q < 7rm, 
the dynamical structure factor S^^{Q,U!) at n = 1 can be expressed as the dynamical 
density-density correlation function of the Sutherland model with coupling parameter 
P = 2 [HI [HI- In order to give this expression, we assume the positive magnetization 
fh without loss of generality, where m = n| — with = N^/N [N„ being the 
number of electrons with cr-spin). In the thermodynamics limit, we have the following 
expression: \ 

5-(Q,c.) = ^ 



/'7r(l — m) /'7r(l— m) 




(4) 

I For the expression of 5^^(Q,t<j) in Ref. 14 , the authors erroneously typed the integration ranges 
of the spinon momenta qi{i = 1, 2) as < (jii < irn^. They should read < g,; < 7r(l — fh), as shown 
in Eq.0. 
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where es{q) is the spinon spectrum: es(<z) = q{vs — q), and ea(q) is the antispinon 
spectrum: e^iq) — q{vs + 5/2), where Vg = 7r(l — m). The purpose of this paper 
is to prove that the analytical expression of S'^^(Q,w) away from half-filling is the 
same as Eq. (0J), if < Q < min[7r?7i, TrnJ. This yields a mathematical proof of the 
strong spin-charge separation in magnetic- field dynamics. We stress that the use of 
the replica type technique is crucial for calculation of the matrix element in S^^{Q, lu). 

This paper is organized as follows. In the next section, we introduce the SU(2,1) 
Sutherland model as an auxiliary. As in the previous study of thermodynamics 0] 
El El Ej and dynamics jSl El El E|i we take the limit /9 — > oo of the coupling 
parameter in order to obtain the analytical knowledge of the l/r^ supersymmetric t-J 
model. The eigenfunctions of the Sutherland model can be expressed in terms of Jack 
polynomials. We discuss the basic features of the Jack polynomials. In sectional we 
derive the matrix element of the dynamical spin structure factor based on the replica 
type technique ^^l- In sectional we present the analytic expression of S^^{Q,oj) 
for finite systems. Section El is devoted to summary. In {Appendix Al we derive the 
dynamical charge structure factor N{Q,uj) in the same method. In |App cndix I 3j we 
show the comparison with numerical results for small size systems [n|. In |Appendix C| 
we present the results on the static structure factors S^^{Q) and N{Q). 



2. Sutherland model with SU(2,1) symmetry 

In this section, we introduce the Sutherland model (201 Ell 1221 ESI with SU(2,1) 
symmetry |24) . and review the basic properties. 



2.1. Notation 

We follow the notations of Refs.|23 EE EZj- For a fixed non-negative integer n, let 
An = {?? = ?72, ■ ■ • , Vn) I Vi ^ Z>o, 1 < i < n} bc the set of all compositions with 
length less than or equal to n. The diagram of a composition ry = (jji, ri2, • ■ ■ , f?n) G A„ 
is defined as the set of points {i,j) € such that 1 < j < rji. The weight ||?7|| of a 
composition ry — (?7i,?y2, ■ • • iVn) € A„ is defined by ||77|| = X)"^! "^he length Z(r/) of 
rj is defined as the number of non-zero rji in rj. The set of all partitions with length less 
than or equal to n is defined by A+ = {A = (Ai, A2, • • • , A„) G A„ | Ai > A2 > • • • > 
A„ > 0}. We also denote a partition A by a'"°6'"''c™'^ • • • or by (a™»,6"Sc'"% • • •) 
with a>5>c>--->0 where is the number of parts which are equal to i. The 
conjugate partition A' of a partition A is a partition whose diagram is the transposition 
of the diagram of A. Hence A^ is the number of nodes in the i-th column of the diagram 
of partition A. In particular we have \\ = '(A). We define the subset A+'^ of the set 
A+ by A+'> = {A = (Al, A2, • • • , A„) e A„ I Al > A2 > • • • > A„ > 0}. Notice that for 
any element A € A+'-^, there exists unique partition /i G A+ such that A = /i + 5{n) 
with 5{n) = (n — 1, n — 2, • • • , 1, 0) S A+'''. For two distinct partitions A, /z G A+, 
we define the dominance order A < /i if ||A|| = ||/i|| and X^iLi ^ XliLi Mi for all 
k = 1, • • • ,n. For a composition 77 e A„, 77+ denotes the unique partition which is a 
rearrangement of the composition r]. Now we define a partial order -< on compositions 
as follows: for iy,r] € A„, v ^ rj ii < t]^ with dominance ordering on partitions or 
if v+ = ?7+ and I]i=i ^« - J2i=i "^i aU fc = 1, • • • , n. 

For a given composition r/ = (771 , 772, ■ ■ • , r]n) and s — {i, j) in the diagram of the 
composition r/, we define the following quantities: 

a,,(s) = ??i - j, (5) 
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a'^{s)^j-l, (6) 
Ivis) = #{fc e {1, • • • , « - <Vk + l<V^} 

+ #{ke{i + l,---,n}\j <7jk<m}, (7) 
=Mke{l,■■■,^-l}\Vk>V^} 

+ #{ke{i + l,---,n}\r]k>m}. (8) 

Here, for a set A, #A denotes the number of elements. The quantities a^(s), a^(s), 
lr){s) and l'jj{s) are cahed arm-length, arm-colength, leg-length, and leg-colength, 
respectively. Since Z^(s) does not depend on j, we also denote it as /^(i). Notice 
that for a partition A S A+, we have 

lx{s) - a;. - z, (9) 

(10) 

Further, for a composition t] G A„ and real parameters r and 7, we define the 
following quantities: 

" (11) 

(12) 

(13) 
(14) 
(15) 



2. 2. Sutherland model and Jack polynomials 

Following R,efs. |14lll7| . we formulate the dynamical spin structure factor S^'^{Q,lu) of 
the l/r^ supersymmetric i- J model based on the freezing technique |15l llfij . 

As an auxiliary, we introduce the Sutherland model |2UI I21[|^l23j with SU(2,1) 
supersymmetry |24| : 

I 0^ 1 /^x2 /3(/3 + A,-) 

^^^= -^2.9^ + I7lzJ Z.,i,2.( )• (16) 

The system has A'h holes, up-spin electrons and down-spin ones, whose 
coordinates are represented by a:^ for i-th hole, x\ for i-th up-spin electron and 
x\ for i-th down-spin electron. We arrange them as x = (xi, 2:2, • • • , xat) = 
[x^, - ■ ■ ,x^^,x\, - ■ ■ ,x\^^,x\, - ■ ■ ,x\^ = {x^.x^.x^). Here the graded exchange 
operator is defined as 

P^,^Y.^fx|-Q,, (17) 

where AC^" is the Hubbard operator which changes from state a to /? one at site j, 
with a,/3 being either h (hole state), or one of a =t, i- The sign factor Qp is —1 if 
/3 =h and 1 otherwise. In order to reproduce the lattice model, we take the limit of 





n(s(^) 


- W^(s) +7), 








d,,(r) = 


n ("''(*) 


+ l + r(/„(s) + l)). 








<(0 = 


n(«''('^) 


+ l + r/„(s)). 












+ r(?,,(s) + l)). 








[0]; = 


n « 


(s)-r/;(s)). 
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large /3 and M, keeping the ratio t — (3/M fixed. Then the particles crystallize with 
equal distance from each other, and the resultant dynamics excluding phonons and 
uniform motion of the center of gravity is mapped to that of the t-J model given by 
Eq.JQ). It can be shown that the intensity of the phonon correlation is smaller than 
the spin correlation by a factor of 0{(3~^). Here we take the lattice parameter L/N 
as the unit of length. Then we have the following relation 

Hcs^tY^D^'^j- (18) 

For fixed numbers of {Ny^, N-^ , N^), the right-hand side of the above relation is the 
t-J model given by Eq.QJ with a trivial constant shift. Note that the symmetry 
of the wavefunction leads to the relation SijPij = —1, where Sij represents the 
exchange operator of the coordinates of particles i and j. The interval / — [l,N] 
denotes {i e Z\l < i < N}. We define 4 = [1, A^h], h = [Nb + 1, A^h + iVJ and 
J| = [A^i^ + iV| + 1 , A^] . The wave function of the ground state for a set of (iVh , A''x , A''^ ) 
is given by 

/ \ /3/2 / \ 1/2 

n n n O-f ■ 

where the complex coordinates z = (2:1, • • • , zjy) are related to the original ones 
X = {xi, ■ ■ ■ , xn) by Zj = exp(27ria;j/L). The spectrum of the Sutherland model 
is conveniently analyzed with the use of a similarity transformation generated by 

o=u (i-f)'''n-r^"^''- (20) 

i^jei ^ ' ' lei 
The transformed Hamiltonian H = O^^HcsO is 

^-^[t) ■ (21) 

Here di is called the Cherednik-Dunkl operator[2Hll2n| and is given by 



di = Zjtt- + /^y^ — - — (1 - ^ik) 
dzi f^. Zi" Zk 



-^^{l-s^k)+P{l-t). (22) 

. , , Zi Zk 

It is known that di can be diagonalized simultaneously by homogeneous polynomials. 
In terms of the monomial z'^ = z'^^ ■ ■ ■ z'^ , the resultant eigenfunctions E^{z] (3) can be 
expressed as Eri{z; (3) = 0''+ lower terms {triangularity)^ and are called nonsymmetric 
Jack polvnomials|3l)l 13 1| . Here "lower terms" means a linear combination of the 
monomial such that v -< ?/. The eigenvalue f]i of Eri{z; P) for di is given by 
m=r]i- Pl'^{i) for z = 1, • • • , A^. 

Since we are dealing with identical particles, the eigenfunction should satisfy the 
following conditions of the SU(2,1) supersymmetry: 

(i) symmetric with respect to the exchange between z^'s; 

(ii) antisymmetric with respect to the exchange between zf s with the same a. 

By taking a linear combination of E\{z\ /3), we can construct a polynomial K\{z\ (3) 
with SU(2,1) supersymmetry 1^ El The above triangular structure of 
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E\{z; (5) is inherited to K\{z\ (3). We specify the set of momenta as A = (A'\ A^, A^) G 
Ajv, where A^ € and A'^ € ^jv'^ (''' =T:i)- For the ground state, we have 

A = Ags^= (A^s,aJ,s'aJ,s) with A^s = (^)^N A^^s = (^(^i) + and 
Acs ~ "^(-^t) (s^^ Figure nj. K\{z;/3) is normahzed so that the coefficient of the 
monomial is unity. 

We define the inner product of functions f{z) and g{z) in n complex variables, 
z = {zi, Z2, ■ ■ ■ , Zn) as follows: 



/ii/|.,|=i 27riz, 



l<fc<i<n 



(23) 



where f[z) denotes the complex conjugation of f{z). We give some examples of 
the SU(2,1) Jack polynomials. For the case of A^h = A^, the polynomial K\(z;f3) 
(A S A~^) reduces to the (monic) Jack polynomial P\{z,j3). The antisymmetric Jack 
polynomial P\{z,l3 + l)nj<je/(^» ~ ^i) S^^^'^ ^y AT^^^^j^^ (z, /3) (A G A^) with 
(A^h, A^i) = (0,0). The SU(2,1) Jack polynomials are orthogonal with respect to the 
above inner product |34l I35| : 

T[NP + 1] N^\N^\N^\ d',{(3) M(3; 1 + (3N) 



{K\, K^)1; = 6\^Z 



(24) 



T[P+1]N p^{f3) dA(/3)/A(/3;l + /3(A^-l))' 
for compositions A, /i G A^ x A^'^ x A^'^ C Ajv- Here px{P) is given by the product 

pxiP) ^ p\mpimpiw with 

A. -\j+f3 



-j-j- Xi- Xj 



(3 



Ai" — \j 



i<y =T,i)- 



(25) 
(26) 



N^ 



2 2 2 

-4 ► < ► 


2 1 
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Figure 1. The diagrams of the ground state with (TV, Ni,,Ni,N^) = (26, 4, 7, 15). 
Acs = (4*, 4^015) + (3*, 5(7), 5(15)) (left) and i^gs = {3^,5(7)) (right). 
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N, - ^'^ A'^ - 1 N j-l .V| - iV| Ni 1 A^i - 1 




Figure 2. The diagrams that appear in the expansions i27V and ^ under 
the conditions lijUl and lijll (the small momentum region). These diagrams 
correspond to the case with {N,Ni„N^,N^) = (26,4,7,15). A = Aqs + 
(i/h,ya,O^T) (left) and fj. = /^gs + (i''',!/") (right) with = (1,0^) and 
u'^ = (1^,0^). The diagrams Aqs and /iqs are shown in Figure [Tl The diagrams 
shown in this figure do not contribute to S^^{Q,lj), owing to = (l,0^h~^) 
(see Eg. 1^1. 

The operators = X°°e'Q^ /VN and = J2j ^j"'^''^' 
expressed for Q = i-Kin/N as Uq — p^/^/N and Uq — p'^/^/N, respectively. Here 
we have introduced power sums — J^iei ('^ ^ ti:?!-!)- In the lattice model 
we have the completeness relation X]<j=| i-^t'^ + -^T' — 1- Therefore, in order to 
calculate S^'\Q^uj)^ we need to know two types of the expansion coefficients, and 
c|, which appear in 

Pm^Aos = E^Ai^^(^;/3)' (27) 
A 

pi-?fAos = E^i^^(^;/3)- (28) 

A 

Using these coefficients, in the lattice limit (/3 cxj), the spin operator 5*1 can be 

expressed by S] = {Xf - Xj^)/2 = f/2 - Xf - Xf. Therefore for Q > 0, we 
have the following relation 

«"w-)4i:('i4)^(£^^'— 

For (5 = 0, there is a finite intensity only at a; = 0, which is given by iVm^/4. It is 
difficult to derive and cj^ for general values of Q = 27:^1/ N . In the next section, we 
show that there occurs a drastic simplification in the small momentum region. 

3. Matrix element 



In the small momentum region, we have some special properties that the calculation for 
the expansion coefficient and norm can be essentially reduced to those of symmetric 



Exact spin dynamics of the SUSY t-J model in a magnetic field 



8 



Jack polynomials 27 . First we summarize the necessary formula to evaluate the 
expansion coefficients c| and c\. Next we obtain the expansion coefficients by use of 
the replica type technique P^. 

3.1. Small momentum region 

We consider the case where the following two conditions are satisfied: 

m < (TVt - ^i)/2, (30) 
m<{N^~ f)/2. (3f) 

(Note that these conditions constitute the small momentum region.) fn this region, 
owing to the triangular structure of the polynomial Kx{z,f3), the composition A 
contributing to the summation (|29|l is restricted to the form A = Acs + ('^'N i^^, O^t) 
(see Figures n and 01 . We define /zgs = (MgS'Mgs) ^'^^^ A*GS = (( ^V^ )^**) ^"^^ 
l^GS ~ ^i^l)- Under the conditions H3()|l and (|3f |l . by extending a calculation by 
Baker and Forrester we obtain the following relation: 

Kxiz;l3) = K^{S;l3') J] zf^ J] (z.-Zj). (32) 

where /?' = /3/(/3 + f), m = (m^s-M^s) + ('^^^') ^ x A+;> and ~z = (z^zi). 
Here K^{z\[3') is a Jack polynomial with SU(l,f) supersymmetry, which is a linear 
combination of A'h + variables non-symmetric Jack polynomials Er,{z\(3'). The 
polynomial -E,,(z; /3') can be obtained by substitution Zj ~ for j S [N\i + Ni + \,N] in 
N variables non-symmetric Jack polynomial E,-i{z;P'). The SU(1,1) Jack polynomial 
Kfj,{z;P') is symmetric with respect to the exchange between z,|^'s and antisymmetric 
with respect to the exchange between zj^^s. For the composition /ios, the SU(1,1) Jack 
polynomial Kf^^g{z, (3') is independent of the parameter and is explicitly given by 

^pos(5,/3')= n^r^'^^' n (^^"^^■)- (33) 

Under the conditions (|30|l and (|31() . the norm of the above states can be reduced to 
the following form ^2 7 : 

_d'AP") fAP"\P"N^) 



/i.H(/3")/.h(/3";l+/3"(iVh-l)) 



hAP') fAP'\ 1 + /3'(A^i + /3"iVh - 1)) 

where (3" = (3'/{(3' + I) = /3/(2/3 + 1) and f3' = P' + I ^ (2(3 + l)/{f3 + 1). We 
remark that for monic symmetric Jack polynomials P\{z,P), the following relation is 
obtained'^: 

(Fa,Pa)^ _ rf',(/3) fxif3;f3N) ^g^^ 



(1,1)^ /iA(/3)/A(/3;l+/3(iV-l))- 

The relations (jSU and (EH) hold if both the conditions < A^ < - 1 for j e /h U 
and A^^ = Xq^ are satisfied. However if Eas.lpUll and lpT|l are not satisfied, the Kx 
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without satisfying Ea. H32|) may involve in the summations (|27f) and (|28|l . In the 
strong coupUng limit (/? oo), the parameter /?' = + 1) approaches unity. 

In this limit, the eigenstates satisfying the above conditions "0 < Aj < — 1 for 
J e /h U Ii and = Ajjg" are SU(2,1) Yangian highest weight states (YHWS) 
in the supersymmetric t-J model, which can be mapped to that of the 

SU(1,1) Sutherland model with coupling parameter /3 = 1 |37|. As employed in 
Refs. [13i38.i.SS«i if the excited states are restricted within the YHWS, the derivation of 
the correlation functions can be reduced to that for the Sutherland model with SU(1,1) 
supersymmetry. We would like to stress that if the small momentum conditions 
l|5n)l and are not satisfied, non-YHWS may contribute in the excited states of 
5-(Q,c.). 



3.2. Replica type technique 

Using the replica type technique |19| , we derive the analytic formula of the coefficients 
and c^. For indices a — h and |, we define the quantity Z"{6) as follows: 

Y[{l-e-'Oz,). (36) 
For any real parameter u, we have the following relation: 



(Z"(6l))" = exp U ^ ln(l - e^^^Zj) = exp -u ^ ^-imePrn 

\ jelc ) \ ra=l ^ ) 

By use of the above relation, the expansion coefficients c" (a = h, j) in Eas. l27|l and 
l)28|l are given by 



, (37) 

0=0 



{Kx,Kx)% u^oiude {Kx,Kx)^ 

where momentum conservation 1 1 A| | — 1 1 Acs II — mis satisfied. We introduce the similar 
type of expansions for the SU(1,1) Jack polynomials: 

pI^K,^, {~z, (3') = ^')' (38) 

/J 

pJ„if^<,3(i,/3')= ^5iif^(5,/3')- (39) 

In a similar manner, the expansion coefficients of c" {a — h, J,) in Eas. H38|l and (|39|l 
are given by 



= lim — 



11 d {Kf.AZ"iO)rK^os)N,+N, 



(40) 



where ||m|| - ||mgs|| = m. 

Under the conditions H3U|) and H31|) . by use of Ea. H32|) . for given compositions 
A — Ags + {i^^, 1^^, O^t) and ^ = ^qs + §, we can show the relation | 

c^ = c^. (a = h,i) (41) 

§ The i/^ and are the partitions with i/^ G A'j^^ and g ^jVj^ ' respectively. 

II In the strong couphng hmit /3 ^ oo, we have 1. For YHWS, Eqs. I34i and 1411 reflect the 

equivalence between the freezing approach and mapping of the eigenstates of the supersymmetric 
J model into those of the SU(1,1) Sutherland model. 
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This relation means that the expansion coefhcients for the SU(2,1) Jack polynomials 
can be expressed by those for the SU(1,1) Jack polynomials, provided the conditions 
H3U|I and (|31|l are satisfied. 

Next, for given parameters {p,q), we consider the following expansion: 

{Z\0)y {zHO)yK,,,iz,P') = ^Xp(/3')^m(^;/3')- (42) 

The formula of the expansion coefficient x^(/3') for arbitrary /i has not been proved 
yet. However, in the small momentum region, we have a formula for XtJ.iP')- For 
A* — f^GS + {i^^, i^^), we have the following relation [77| : 

By use of Eas. H4()|l and l|43|l with (p, q) — (0, u), we can obtain the coefficient c|. For 
A = Acs + (^^^ O^T), it is given by 

where to = + jji^'^ll. In a similar manner, we can derive the coefficient (see 

Ea. ljA.3|) ). The first term in Ea. H44|) can then be represented by — /3"c^, which vanishes 
unless = (0^^). We consider the quantity c\ = c\^+ /^''^aj which is nothing but the 
second term of the right-hand side in Ea. l(^ . Owing to the factor f^h{P"; 0), the c\ 
vanishes unless A = (0^''). Therefore, the cx can be simplified as 

^A^IKII^T^x W~H)- (45) 

We notice that while is related to the expansion coefficients for the SU(1,1) Jack 
polynomials, the quantity c\ given by Ea. (|45|l is related to those for the monic 
symmetric Jack polynomials. In fact, for the monic symmetric Jack polynomials, 
one has the following formula |4()j : 



In the strong coupling limit (/3 — > oo), we have /3" 1/2, and therefore the c\ 
approaches + cJj/2. This is precisely the quantity which appears in S^^{Q,A 
(see Ea. (|29ll '). Thus, no charge excitation contributes to S^'^{Q,A in the small 
momentum region. Since we have (3' 2 in the limit P oo, the coefficient cx 
vanishes owing to the factor [0]^a, if partition ly^ contains s — {i,j) — (2,3). Then, 
S^^{Q,A is determined by three parameters (Aa, Agi, As2). Partition i^^ is restricted to 
u"" = (Aa,2^=2~\l'*'=i"^=^0^i~-^=i). These three parameters are related directly with 
momenta of the elementary excitations: one antispinon and two spinous. 



4. Results 



Using the relations (|34|l and (|45|l . we can express the S^'^{Q,A in terms of three 
parameters (Aa, Asi, As2). In the small momentum region (see Eas. (|30ll and lj2U), we 
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Figure 3. 5^^(Q,aj) for {N,Ni„Nf,Ni) = (16,0,10,6) (left) and (16,2,9,5) 
(right). Each spectral weight is proportional to the area of the oval. For dark 
shaded ovals in the right figure, excitation energies and spectral weights agree 
with those in the Haldane-Shastry model (the left figure), which can be expressed 
by Eq.OZj and l48l . 




Figure 4. S^^(Q,u}) for (N, Ni,, , N^) = (16,0,12,4) (left) and (16,2,11,3) 
(right). For dark shaded ovals in the right figure, excitation energies and spectral 
weights agree with those in the Haldane-Shastry model (the left figure), which 
can be expressed by Eq. I4YI and I4SI . 



obtain the following expression: 



2N 



A3i>A,2,Aa 

X (Aa - l)(Aa + N-{Nt- N^) - l)(A,i - A,2 + \) 



n 



L (Aa + 2Asj -J- l)(Aa + 2Asj - ]) 



j = l -L [■^sj 2~J^ It 2 -^sj + —\ 



(47) 



where Q = 2'Km/N and = Aa + '^j^iiKj — !)• The excitation energy AE\ is 
given by 

2 



AEx/t 



(Aa 



N Nt - Ni Aa 
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iV iVt - iVi , 1 



(48) 



In the case of As2 — 0, where only one spinon is excited, we need to modify the above 
results as follows: 

. 1 r[i]r[A.i + i]r[f - ^ + i]r[f -^^x., + i] 
^^'"^ 2^N r[A.i + i]r[f - ^ + i]r[f ~ ^ - + 1] 

X S{uj - AEx), (49) 

where Q — 2it\si/N. The excitation energy /S.E\ is given by (Asi, As2, Aa) = (Asi, 0, 2) 
in Ea. H48l) . We have checked the validity by comparison with numerical results up 
to iV = 16 [12 (see [Appendix B) ) . In Figures [3| and 0] we show the results for 
A'^ = 16. From comparison with numerical results, the analytic expression of the two- 
spinon plus one-antispinon contribution can be applied in the wider range of {Q,uj) 
(see [Appendix B[ ). The analytic expressions H47|l to 1)491) coincide with those for the 
Haldane-Shastry model From the above results for finite systems, we can derive 
the analytic expression of S^'^{Q,ui) in the thermodynamic limit (see Eq.lQ). Note 
that the contribution in the case of As2 = vanishes in the thermodynamic limit. Thus 
we have proved analytically that in the momentum region < Q < min[7rm, fcpj], the 
structure factor S^''{Q,uj) is not affected by hole doping. Here /cf,o- is given by Trflo-. 
In this region, S^^{Q,u;) diverges as (uj — CsiQ))'^^"^, as the frequency approaches the 
lower edge corresponding to the spinon dispersion uj — es{Q). The obtained S^^{Q, uj) 
has the same form as the dynamical density-density correlation function of the spinless 
Sutherland model with coupling (3 — 2 except the momentum range [^ [32j- For 
< Q < min[7rm, fep^x], the static structure factor S^^{Q) can be evaluated by the 
integration J dLuS^^{Q,uj), which reproduces the expressions presented in [Appendix C[ 



5. Summary 

By use of the freezing technique on the Sutherland model with SU(2,1) supersymmetry 
and the replica type technique, we have obtained the dynamical spin structure 
factor S^^{Q,uj) with Q < min[7rm, A:f,x] in the supersymmetric t-J model with l/r^ 
interaction. The 5^^(Q,a-') has the same form as that of the Haldane-Shastry model 
in this small momentum region. In the thermodynamic limit, two spinous and one 
antispinon contribute to 5^^((3, oj). Therefore, S'^^'{Q, us) is not affected by hole doping 
in this region. Thus we have proved the strong spin-charge separation in S^^{Q,uj), 
which was numerically obtained in the previous paper |12| . From comparison with 
numerical results, we have found that the analytic expression of two-spinon plus one- 
antispinon contribution can be applied to the wider range of (Q,a;). 
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Appendix A. Dynamical charge structure factor N{Q,uj) 

We derive the dynamical ciiarge structure factor N{Q, uj) for Q < kF,i H by use of the 
replica type technique JHl- The N{Q,uj) is defined by 

7V(Q,c^)= J2\{,y\nQ\0)\^diLO-E, + Eo), (A.l) 

where ng = ^Ylii'^i^^'^^ I ■ ^.s in our previous studvfTT]. using the expansion 
coefficient c'jj in Ea. (P7|l . the dynamical charge structure factor N{Q,lu) can be 
expressed as 

N[Q. -) = ^ E ('^a)' 7#^#V'5(- - ^E,). (A.2) 

for the momentum Q > 0. Using the replica type technique '19', we can obtain the 
coefficient c^. By use of Eas. l4()|l and 143|l with (p, q) = {u, 0), the expansion coefficient 
can be derived as 

in the small momentum region m < (Ni — l)/2. Here we have used the property that 
/y=(/5';0) vanishes unless = (0^^). In fact, as shown in Ref.|12|, the excited states 
contributing to N(Q,lu) in this small momentum region are restricted to the case 
where A = Acs + (j^^, O^i, O^t). in this case, the SU(2,1) Jack polynomial Kx{z,l3) 
in has a form Kx{z,(3) = P^h(z^,/3") x Kx^^. Therefore the coefficient 

can be derived via Ea. H46|l as well. Namely, the use of the replica type technique is 
not essential for the derivation of in contrast to the case for c\^. The norm can 
be evaluated by the reduced formula Ea. (|M|l . In the strong coupling limit /3 ^ oo, 
the coupling parameter /3" becomes 1/2. In the thermodynamic limit, the dynamical 
charge structure factor N{Q,uj) for < Q < fcp,! can be expressed as [T7] 

X 6 [ Q - gall - y ]qj \ 5 \ UJ ~ eah(gah) - f^h{qj 







dqaii / 


dqi / dq2 




Jo 


2 N 




- gah - E ij 


\ 6 \uJ~ eah(gah 


3 = 1 J 





\qi - '72|eah('7ali) 



11^=1 Vehfe) (2gj + gah)^ 



(A.4) 



'■3 = 

where the Fermi momentum fcp is given by fcp — 7rf;,/2, eh{q) is the holon spectrum: 
^h{q) — q{vc + q) and ea.\i{q) is the antiholon spectrum: eah(<z) — q{vc — <z/2). Here 
the charge velocity Vc is Vc = 7r(l — n). This expression has the same form as the 
dynamical density-density correlation function of the spinless Sutherland model with 
coupling parameter (3 — 1/2 . 

Appendix B. Comparison with numerical results 

We make a comparison between analytic results and numerical ones in S^^{Q,uj) 
[H) . In Tables IbTI and lB2l we present the cases (iV, iVh, A^f ^i) = (16,2,9,5) 

^ For the derivation of N[Q,u)), the condition Q < 7rm, i.e., Ea. l3Ui is not necessary. 
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Table Bl. Comparison between analytic results and numerical ones 1121 for 

{N,Ni„N^,Ni) = (16,2,9,5). 



(Asi,As2,Aa) Q/tt u/t I\ (analytic) I\ (numeric) 

(1,0,1) I ^(~ 0.848169) ^(~ 0.03409090) 0.03409080 

(2,0,1) I ^(~ 1.38791) ^(~ 0.0505050) 0.05050442 

(1,1,2) ^(- 1.85055) ^(- 0.02462121) 0.02462107 

(3,0,1) f ^(- 1.61923) ^(- 0.06926406) 0.06926359 

(2.1.2) ^(- 2.39029) ^(- 0.036931818) 0.03693162 

(1.1.3) ^(- 3.00715) 3|2(- 0.01988636) 0.01988634 

(2,2,2) %-(- 3.08425) gi|o(- 0.00200617) 0.00200648 



Table B2. Comparison between analytic results and numerical ones 1121 for 
(N,Ni„NT,Ni) = (16,2,11,3). 



(Asi, As2, Aq) 


Q/tt 




Uj/t 


h 


(analytic) 


I\ (numeric) 


(1,0,1) 
(2,0,1) 
(1,1,2) 


1 

8 
1 
4 

3 


128 V- 
Stt 

64,^- 
8 V- 


t 0.539744) 
! 0.771063) 
- 1.2337) 


28 V— 
35 V— 
112 V- 


: 0.03571428) 
0.057142857) 
i 0.02678571) 


0.035714172 
0.057142144 
0.026785570 


(1,1,3) 


8 
1 


277r^ (, 
128 ^" 


- 2.08187) 


224 V— 


0.022321428) 


0.022321422 


(1,1,4) 


2 

5 


57r^ 
16 V- 


- 3.08425) 


560 V— 


0.019642857) 


0.019642881 


(1,1,5) 


8 


557r^ (, 
128 V 


- 4.24085) 


56 (— 


: 0.01785714) 


0.017857178 



and (16,2,11,3), respectively. Our analytic proof is restricted to the case where 
Q < min[fcF_|, ttto]. However, the analytical expression of the two-spinon plus one- 
antispinon contribution can be applied in the wider range. As a result of hole doping, 
the integration ranges of the spinon momenta in Eq.Q) are changed to < < /cp.i 
for i = 1 and 2. From comparison with numerical results (rnill2) . we find that a similar 
fact occurs also in the N{Q,uj). Namely, although analytic derivation of N{Q,uj) is 
restricted to the region < Q < kp^i, the expression of the (right-moving) two-holon 
plus one-antiholon contribution can be extended to the integration range shown in 
Eq.E3I). 

Appendix C. Static structure factors 

We consider the static structure factors S^^{Q) and N{Q). There are several ways 
to obtain these quantities. If one knows the dynamical structure factors S^^{Q,lS) 
and N{Q,iij), then the static structure factors can be obtained by integration over 
uj. Gebhard and VoUhardt calculated the static structure factors for ?7i = 0, from 
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the Gutzwiller wave function ,.43^ . For general ?7i, Forrester derived the analytic 
expressions of the equal-time two-point correlation functions ^1 021 ^- In the 
following, we obtain the the static structure factors by Fourier transformation of these 
equal-time two-point correlation functions. 

The equal-time two-point correlation functions are defined by C^^{x) = 
(0|S'^5^|0) and C'^^ix) = {O\n^no\0). They are expressed as follows, 



- 9 


n — m? 


T ' 


4 




1 

_l_ 






4 





dx 



+ (1 - 4.o) 



dx 



Sc{x) 



dus^{u) 



dus^iu) 



where s_(x) is s-.{x) = Ss{x) — Sc{x). Sa{x) (a —c and s) are given by 
, , sinVaX 

Sa[X) = . 



(C.l) 



(C.2) 



(C.3) 



By Fourier transformation we obtain the analytic expressions of S^^{Q) and N{Q). 
Taking into account of the Umklapp process, we obtain jJH] 

5-(Q) ^ _± ^ + SiiQ) + 5/(27r -Q) + 5//(Q) + 5//(27r - Q), 

47r 

N{Q) = - + NiiQ) + Nj{2n - Q) + Njj{Q) + Njj{2n - Q), 

for momentum < Q < 27r. The Si{Q) is given by 

1, 



SiiQ) 



dx cos Qx 



[Ss{x)Y 



9{2v, - Q) 



Q_2vs 
An 



■In 



dx 
Q 



Ss{x) 



duSs{u) 

(C.4) 



where 0{x) is 9{x) = 1 for positive x, and otherwise. This contribution has 
the same form as the level-level correlation of the random matrices for symplectic 
ensembles |23l 1461 14 7| . In fact, the S^^{Q) of the Haldane-Shastry model can be 
expressed by Vs/{2Tr) + Si{Q) + Si{2n ~ Q) The Sii{Q) contributes for finite 

hole doping (h < 1), which is given by 



SiiiQ) 







dxcosQx —Ss(x) 
' dx 



Vc 

47r' 



• + Vc 



— In 

Stt 



0, 



dusc{u) 

for < Q < Vs ~ Vc, 

, for Vs - Vc < Q < Vs 
for Q > Vs + Vc- 



The divergence at Q = 7r(l — m) in Si{Q) is removed by hole doping. The static 
spin structure factor has the same form for the Haldane-Shastry model in the region 

+ In Ref. l44l . the normalization factors are different. When the spin correlation is divided into 
C^^(a;) = C^^{x)/'i + C^^{x) + C^^{x), the correlation functions C^^{x), C^^{x) and C^^{x) 
correspond to p'^hoo{x), PoPih^iix) and p^^h^^{x) in Ref. 1441 . respectively. Tractable expressions 
for the correlation functions were derived by Kuramoto l45l . 
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Q < 2kp i. This region contains "Q < min[7rm, fcp.j,]" , where the S'^^{Q, lu) of the l/r^ 
supersymmetric t-J model has the same form as that of the Haldane-Shastry model. 
For the momentum < Q < min[7rTO, fcp J, the w-integration of the S^^{Q,uj) (see 
Eq.l^J) reproduces the above expression. 

Next we consider the static charge structure factor N{Q). The Ni{Q) is given by 



Nj{Q) = 2 [ dxcosQx ~[scix)]'^ — ( -i-Scix)] [ dusc{u) 
Jo I \dx ) 



+ — In 



Vc Q , 

m 

TT 27r 



Vr 



Q + Vc 



, for < Q < 2vc, 
for 2vr < Q. 



(C.6) 



Q - Vc 

This term has the same form as the level-level correlation of the random matrices for 
orthogonal ensembles |46l E7|. In fact, the static structure factor of the Sutherland 
model with coupling parameter /3 = 1/2 is given by Uc/tt + Ni{Q) '23'. The Njj{Q) 
is given by 



Nii(Q) = 2 dxcosQx [ -r^cix)] / dus. 
Jo \dx J J^ 



0, 

Q , Q + Vc Q -Vs- 

— m — 

2tt zvr 



for < Q < Us — ^'c , 
, for Vs - Vc < Q < Vs + V, 

„ — „ , for Q>Vs + Vc. 

2tT Q - Vc 

For < Q < kp^i, the w-integration of the N{Q,uj) (see Ea. (IA.4|) ') reproduces the 
above expression. 



Vc 
TT 



Q , Q + Vc 

■ m ■ 




S 1/2 



n 




Figure CI. The shaded region shows the permissible one (0 < n < 1 and 
< < n) of the (n,fn) space in the t-J modeh This region can be classified 
into six regions with different expressions of S'^^(Q) (left) and N(Q) (right). 

We can rewrite S'^^{Q) and N{Q) more explicitly. Because of the reflection 
symmetry against Q = tt, it is enough to consider Q < n. For convenience, we 
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define the following functions: 

X X 



Si{x) 

S2a{x) 
S2b{x) 
Sza{x) 

Siaix) 
Sibix) 



4 

n - 



■In 



1 — 771 



1 — m 
-?ln. 



777 X 

m 

2 8 



1 



777 



77—1 
77 + 3777 



X X 

4-8'° 



I — X 



1 + 777 — a: 



1 



X 



X . 



77 — 2777 + 1 

4 '■ 

77 + 2777 — 1 



In 



■In 



m 
rfi - 



1 - 77 



1 



1 + 777 — a; 



In 



■In 



1-777 



1—777 



54c(x) =---ln- 



(C.8) 



The expressions of S^-^-{Q) can classified into six cases (a) to (f) (see Figure lCTl flcft)'). 
(a) When m > 1/2, we have 



Si {Q/tt) , 

S2a (Q/tt) : 
Ssa (Q/tt) ; 



777j, 
- 777) 



for < Q < 7r(77 — 777), 
for 7r(77 — 777,) < Q < 7r(2 — 77 - 
for 7r(2 - 77 - m) < Q < 27r(l 
for 27r(l - m) < Q < TT. 

(b) When 777 < 1/2 and 777 > —77/3 + 2/3, we have 

S'i(g/7r), for < Q < 7r(77- 777), 

S2a (Q/t^) I for 7r(77 — 777.) < Q < 7r(2 — 77 — 777) 

Ssa {Q /t^) , for 7r(2 — 77 — 777) < Q < 27r777, 
5*46 (Q/tt) , for 27r777 < Q < tt. 

(c) When r?7 < —77/3 + 2/3, ttt. > — n + 1 and ttt. > 77/8, we have 

' SiiQ/ir), for < g < 7r(f7- m), 

S2a (Q/tt) , for 7r(77 — 777.) < Q < 27r7fl, 

Ssb (Q/i") , for 27r777 < Q < 7r(2 — 77 — m), 
5*46 (Q/tt) , for 7r(2 — n — m) < Q < tt. 

(d) When 777 < —77 + 1 and ttt. > f7/3, we have 

' Si{Q/n), for < Q < 7r(T7- T77), 

S2a (Q/ti") , for 7r(77 - ttt) < Q < 27r?7T, 

Ssb (Q/tt) , for 27rT77 < Q < 7r(f7 + ttt), 

S4C (Q/tt) , for 7r(f7 + ttt) < Q < tt. 

(e) When ttt < —77 + 1 and ttt < tt/3, we have 

Si{Q/t:), for < g < 27rTTT, 

S2b {Q/t^) , for 27rT7T < Q < Tr{n — fh), 

Ssb (Q/t^) , for 7r(fT — ttt) < Q < 7r(77 + ttt). 

Sic {Q/t^) , for 7r(f7 + 777) < Q < TT. 



(C.9) 



(C.IO) 



S^^{Q) 



(C.ll) 



S^^{Q) 



(C.12) 



S''{Q) 



(C.13) 
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SiiQ/ir), forO<Q<27rTO, 

5*26 (Q/tt) , for 27rm < Q < n{h ~ m), 

5*36 (Q/tt) , for 7r(fT, — fh) < Q < 7r(2 — h — rh) 

Sib iQ/'^) ; for ■"'(2 — n — fh) < Q < tt. 
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(C.14) 



To describe N{Q) as weh, we define the following functions: 



Ni{x) 

N2a{x) 
N2bix) 

Nsaix) 
N3b{x) 
Niaix) 

Nib{x) 
Nicix) 



■In 



2 

fh 



1-h 



1 



n 



X .T , 1 
- + - In - 
2 2 1 



2n + I In 

X 



n 



m 
1 + a; 



1 — n + .T 
1 + n - 



In 



1 



In 



m 



2 - -n 
2 



a; 



= 2h 

= 2 
= 1- 



m a; 
Y ~ 2 2 

1 — n 



In 



1 — TO 



In 



1 



TO 



2n- 



■In 



In 



N{Q) 



for < Q < 7r(fi — to), 
for 7r(fi — fh) < Q < TT{h + to), 
for 7r(fi + to) < Q < 2-Kh, 
for 27rn < Q < tt. 



N{Q) 



The expressions of N{Q) can be classified into six cases (a) to (f) (see Fi] 
(right)). 

(a) When n < 1/2, we have 

^1 (Q/tt) , 

N2a (Q/tt) 
N3a (Q/tt) 
N4a iQ/n) 

(b) When n > 1/2 and fh < — 3n + 2, we have 

A^i(Q/7r), for < Q < 7r(n-?fi), 

N2a (Q/tt) , for Tr{h ~ fh) < Q < 7r(n + to), 

N3a{Q/Tr), for7r(n + TO)<g<27r(l-n), 

A^46(Q/7r), for 27r(l - n) < g < TT. 

(c) When to > 3n — 2, to < —h + 1 and fh > — 3n + 2, we have 

A^i(Q/7r), for < Q < 7r(n-TO), 

N2a{Q/Tr), for7r(n-?7i) < g< 27r(l-n), 

^36(g/7i-), for 27r(l-n) < Q< 7r(n + ?fi), 

Nib iQ/Tr), for 7r(n + to) < g < TT. 

(d) When to > —h + 1 and to > Bn — 2, we have 

A^i(g/7r), for < g < 7r(n-TO), 

N2a{Q/TT), for7r(n-TO)<g<27r(l-n), 

NsbiQ/ir), for 27r(l - n) < g < 7r(2 - n 

Nic iQ/Tr), for 7r(2 - n- to) < g < TT. 



(C.15) 
rure ICll 



(C.16) 



(C.17) 



A^(0) 



(C.18) 



NiQ) 



to). 



(C.19) 
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(e) When m > —n+1 and 



m < 3ri — 2, 



we have 




m < 3n — 2, 



Ni (Q/tt) , 
N2b {Q/tt) , 
Nsb iQ/n) , 
Nic (Q/tt) , 



we have 



for < Q < 27r(l - n), 
for 27r(l — n) < Q < 7r(n — m), 
for 7r(f;, — m) < Q < 7r(2 — n — rh 
for 7r(2 — n — m) < Q < tt. 



(C.20) 



A^i (Q/tt) , 
A^26 (Q/tt) , 

Nib (Q/t:) , 



for 
for 
for 
for 



< Q < 27r(l - n), 
27r(l - n) < Q < 7r(n- m), 
7r(fi — m) < Q < 7r(n + m), 
7r(fi + m) < Q < TT. 



(C.21) 



In the limit m ^ 0, the above expressions of S^^{Q) and N{Q) reproduce the results 
in Ref.|121. 
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